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Influence of Impeller Shroud 
Forces on Turbopump Rotor 
Dynamics 


The shrouded-impeller leakage path forces calculated 

analvred to answer two Questions. First, because of certain characteristics oj 
results of Childs the forces could not be modeled with traditional approac 
therefore an approach has been devised to include the forces " convent, w 
rotordvnamic analyses. The forces were approximated by tradmonal V t ^ , o ’ d ^ ir f r 
ine and inertia coefficients with the addition of whirl-frequency-depende 
and coupled stiffness terms. The forces were found to be well-modeled with 
this approach. Finally, the effect these forces had on a ^mp/e roror-^ar;ny>5 
was analyzed, and, therefore, they, in addition to seal forces 
Jeffcott rotor The traditional methods of dynamic system analysts were 

to incorporate the impeller forces and yielded results for the eigenproblem.frequenc 

tZ f^Lncyoffree vibration as well as system stability. All results lead to the 
forces have little influence on ^ fequ^ but can hate 
appreciable effects on system stability. Specifically, at hig ' W expected 
at the leakaee Doth entrance , relative stability is reduced. T ^ , 

response characteristics that occurred are attributed to the nonlinearity of the model. 


Introduction 

The vibration of centrifugal pumps has received increasing 
attention recently because of the inability of current analysis 
techniques to adequately predict the dynamic characteristics 
of pump designs. Failure to accurately predict vibrations has 
resulted in the loss of considerable amounts of money in down 
time from severe vibration problems. Massey (1985), for ex- 
ample, described an eleven-stage pump that became unstable 
when its running speed exceeded its critical speed by 25 percen . 
In other words, it whirled at 80 percent of running speed. 
Another example occurred in the High Pressure Oxygen Tur- 
bopump (HPOTP) of the Space Shuttle Main En 6ine (SSME) 
which also whirled at 80 percent of running speed (Childs and 

Moyer, 1985). _ . , 

The hydrodynamic forces generated in many of the flu 
filled gaps within the pump are well-established contributors 
to the problems cited above. This article is concerned specif- 
ically with the forces developed along the leakage path between 
the impeller and the shroud, as shown in Fig. 1. As the fluid 
is discharged from the impeller, some will return to the lower- 
pressure, suction side by way of this leakage path. Leakage is 
minimized typically with wear ring seals, as shown. Hydro- 
dynamic forces are developed along the entire leakage path, 
i e., the shroud section as well as the seal. Note that these are 
only parts of the total force on the impeller and that the 
impeller/volute region, balance drums, inducers, shaft, etc. 
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also contribute to the dynamics of pumps, although they are 

n °Unfortun^lyf measurements of the *ota‘ fora 

SmTdif^ However, some 

StT have been reported that are of importance here. For 
example, research at Cal Tech (Adkins, 1976) led to the con- 
clusion that the pressures in the shroud annulus contribue 
from about 50 percent to 75 percent of the total stiffness acting 
on the impeller. The pump used included separation rings and 
an enlarged shroud clearance space to minimize leakage flo 

f °Bolleter et al. (1985) used a rocking-arm mechanism to ver- 
tically translate the spinning impeller. Their pump had formal 
clearances in the leakage path; and, therefore, the forces meas 
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ured were more realistic. The results demonstrated that the 
nondimensionalized cross-coupled stiffness value was about 
twice that measured on the pump at Cal Tech (Jery et al. f 
1984). This suggests that the leakage path region in pumps can 
also reduce the stability of impeller motion. Bolleter et al. 
(1989) have recently presented extensive force-coefficient re- 
sults at off-design conditions. 

Almost all of the analytical attempts at predicting leakage 
path forces have been concerned with seals. However, Childs 
has extended his previous techniques in seal analysis (1980, 81, 
82a, 82b, 82c) to apply them to the flow within the clearance 
space surrounding the shroud (1987). His work will be used 
to model the seal and shroud forces in this article and will now 
be summarized. 

Childs used a bulk-flow approach to obtain the governing 
equations of the flow of a differential element of fluid. In the 
seal problem, three equations were required: axial- and cir- 
cumferential-momentum equations as well as the continuity 
equation. In the shroud problem, the axial-momentum equa- 
tion was replaced by a path-momentum equation, introducing 
additional terms which described the centrifugal and Coriolis 
accelerations of the fluid element. After a perturbation ex- 
pansion of the equations in the eccentricity ratio, the resulting 
relationships were solved at various values of the whirl fre- 
quency, fi. Integration of the resulting pressure distribution 
acting on the shaft or shroud yielded the radial and tangential 
force as functions of whirl frequency. In addition, solution 
was carried out at various fluid circumferential velocities at 
the shroud entrance. This variable is of primary importance 
in system stability (cross-coupled stiffness). Figure 2 shows the 
results of the shroud problem for three different values of the 
inlet swirl velocity. The leakage path analyzed corresponded 
to that used by Bolleter et al. (1985) which ran at 2000 rpm. 
Figure 3 provides the dimensions of the pump and other data 
of importance. The curves in Fig. 2 define force coefficients 
corresponding to a circular orbit of amplitude and whirl fre- 
quency ratio /. The impeller force components are linearly 
proportional to the orbit amplitude A. The sharp deviations 



Whirl Fre qu ency Ratio (ND) 

Fig. 2(4) Radial force coefficient curves (Childs, 1987) 



Whirl Frequency fUoo (ND) 

Fig. 2(b) Tangential force coefficient curves (Childs, 1987) 


in the radial and tangential force coefficients of Fig. 2 are 
caused by excitation of a “centrifugal-acceleration mode” nat- 
ural frequency. Childs (1990) has presented results for damped 
natural frequencies and eigenvalues for these modes. 


Nomenclature 
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cross-coupled stiff- 
ness coefficient 
which is a function 
of the whirl fre- 
quency ratio ( F/L ) 
length of seal or im- 
peller-shroud path 
(L) 
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direct mass coeffi- 
cient (Childs, 1987) 
rotor whirl radius 
(L) 

clearance at leakage 
path entrance (L) 

clearance ratio 

rotor whirl ampli- 
tude (£) 
inlet radius of 
shroud leakage path 
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nondimensionalized 
tangential fluid ve- 
locity at leakage 
path or seal entrance 
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(inlet swirl ratio) 
(Childs, 1987) 

= rectangular coordi- 
nates of rotor center 
(L) — 

= complex rotor whirl 
amplitude (£) 

= complex rotor coor- 
dinate (L) 

= phase angle of whirl 
from harmonic exci- 
tation 

= rotor whirl fre- 
quency (l/7> 

= natural frequency of 
free vibration (1 /T) 
= shaft speed (1/7*) 

= shroud leakage path 
coefficient subscript 
= rotor coefficient 
subscript 

= seal coefficient sub- 
script 

= shaft coefficient sub- 
script 

= overall coefficient 
subscript 
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Pump da La 

shaft speed = co 5 = 2000 rpm 
inlet clearance * Q 0 = 3.55 mm 
impeller exit angle * 22.5 0 
flowrate = 130 1/sec 
head developed = 68 m 
seven- bladcd impeller 
working fluid is water at 80 °F 

Fig. 3 Dimensions and data for the pump tested by Bolleter et al (Bol- 
leter et al, 1965) 


In the articles cited above, Childs initially suggested the 
following conventional, linear relationship to model the forces. 
It applies for small motion of the rotor about the centered 
position. 



( 1 ) 


If circular whirl at frequency, fl, and radius, Q t is assumed 
such that 

jc=QcosQ/, (2) 

y=Qs\nQt t 

the radial and tangential force components can be stated 

After nondimensionalizing the terms as defined in Childs 
(1987), Eqs. (3) become 

-=Mf 2 -cf-R 

4 (4) 

5--cr+r, 

Q 

where the definitions of the coefficients are included in the 
Nomenclature and the other variables are defined as follows. 


p- nondimensionalized force 

1R l LAP 

d=— % nondimensionalized clearance 

Go 



whirl frequency ratio 
force coefficient. 


The quadratic relationships of equations (4) are to be used to 
obtain the dynamic coefficients of equations (1). Specifically, 
a least-squares curve fit of the force curves will yield the coef- 
ficients. This is an adequate approach for seal forces because 
they, in general, do follow a quadratic trend. However, the 
curves of Fig. 2 illustrate that shroud forces cannot be ap- 
proximated by equations (4); and, therefore, a traditional lin- 
ear model as expressed in equations (1) will not adequately 
represent the forces. Although these results are unexpected, 
experimental measurements made by Franz and Arndt (1986) 
of impellers with inducers yielded very similar results. 


Problem Definition 

At this point, it is not known, assuming the force curves of 
Fig. 2 are valid, how they can best be included in traditional 
rotordynamic calculations of stability, response, critical speeds, 
etc. To solve this problem, the present research has been con- 
ducted to answer the following questions. 

(/) How can the force curves found by Childs (1987) be 
incorporated in the traditional rotordynamic analyses 
of system response, stability, critical speeds, etc.? 

(2) What effects do the forces have on a rotor-bearing 
system with respect to stability, natural frequencies, 
imbalance response, etc.? 

Before answering these questions, observe that the forces pre- 
sent no new problems in a transient analysis because they could 
be applied as whirl-frequency-dependent forcing functions to 
a simple rotor model such as that developed by Jeffcott (1919). 
Although the force components are linearly proportional to 
the orbit amplitude, the equations of motion would be non- 
linear because the forces depend on the whirl frequency, de- 
fined kinematically as 


H xy-yx 
w,~ w ! (x 2 +y z )' 


(5) 


but this could be handled by integrating the equations of mo- 
tion numerically with traditional methods. However, a tran- 
sient analysis alone is not sufficient to obtain an understanding 
of the forces, and this is why answers to the questions posed 
above are necessary. The following section describes the model 
used to analyze the shroud forces, and results are then given 
for the eigenvalue problem and frequency response, respec- 
tively. 


Rotordynamic Model 

As shown in Fig. 4, the model is a Jeffcott rotor (Jeffcott, 
1919) under the application of seal and shroud forces. The 
rotor itself represents a double-suction impeller, yielding two 
symmetrical leakage paths. In summary, the forces on the rotor 
are the shaft stiffness, two identical seal forces, and two iden- 
tical shroud forces. The mathematical representation of each 
force will now be discussed. 

A 90.7 kg (2001b) rotor is assumed. In addition, the shaft 
stiffness used (K sh = .2987 MN/m) makes the natural frequency 
of forward whirl (to be defined later) equal to 80 percent of 
shaft speed. This choice is not completely arbitrary since this 
frequency corresponds to that observed in the examples cited 
in the Introduction (Massey, 1985; Childs and Moyer, 1985). 

The seals modeled are smooth, wear ring seals with a clear- 
ance and diameter of .36 mm and 236 mm, respectively. The 
least-squares approach of Childs described earlier has been 
carried out, and the resulting dynamic coefficients are given 
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seal forces 



shroud forces 

Fig, 4 Jeffcott-based, double suction Impeller leakage path model 


Table 1 Coefficients for one wear ring seal 


UtL 

M, 

kg 

c. 

Ns/m 

Ns/m 

K, 

MN/m 

k, 

MN/m 

.5 

.3111 

3356.0 

78.6 

.6121 

.4628 

.6 

.3398 

2167.0 

89.6 

.5443 

.4611 

.7 

.3714 

970.9 

101.3 

.4755 

4455 


in Table 1. The analysis used has at times under-predicted the 
true value of the direct stiffness. Therefore, K s is doubled when 
the coefficients of the entire model are assembled. This mod- 
ification of K s has no influence on the results, since the shaft 
stiffness is selected to yield a desired system natural frequency. 

As described in the introduction, conventional stiffness, 
damping and mass coefficients cannot completely describe the 
impeller-shroud forces calculated by Childs and shown in Fig. 
2. Therefore, the following method has been devised. First, a 
curve fit of the force curves is carried out, yielding dynamic 
coefficients which describe the quadratic variation of tfieforces 
with respect to the whirl frequency ratio as described above. 
The difference between these curve fits and the actual force 
curves is modeled by direct and cross-coupled stiffness coef- 
ficients which are functions of whirl frequency ratio. The math- 
ematical equivalent of this approach will now be discussed. 

The following modified forms of equations (4) express the 
method with which the impeller-shroud forces are modeled. 



FlO- 5« Dependency ot K’ on whlri frequency ratio 



Table 2 Approximating coefficients for one leakage path 


Mi 

c L 

CJ 1 

K l 

1 ^1 

kg 

Ns/m 

Ns/m 

MN/m 

MN/m 

3.469 

2037 

969.9 

r. 06781 

.3137 


F, 


--MtJ 1 - ? L f- + </)} 


F, 

— w 

<7 




( 6 ) 


where R’ and k' represent the nondimensionaiized deviatior 
between the force curves and the approximating quadrati 
expressions. Figure 5 illustrates R' and k\ respectively. Not 
that, although they represent forces proportional to displace 
ment, they are not traditional stiffness coefficients since the 
are functions of whirl frequency ratio. Separating the fore 
coefficients into whirl-frequency-independent and whirl-fre 
quency-dependent terms does not introduce any approximatioi 
to the force definitions. Equations (1) are now rewritten as 


r 

l-k L ~ 


K*(f) k L + k*(f) 

**(/) Kl+K'V) 


0 - 


In a mathematical sense, the tangential force accounted for 
by k * could be modeled as a damping coefficient; however, 
the result of Fig. 5(b) confirms the wisdom of modeling the 
deviation of the tangential-force coefficient by a whirl-fre- 
quency-dependent cross-coupled-stiffness-coefficient. Specif- 
ically, k * oscillates about zero and approaches zero as the inlet 
swirl ratio is reduced. Finally, observe that equations (7) are. 


in general, nonlinear differential equations. To understand this 
statement, recall that the general definition of the whirl fre- 
quency ratio,/, is given in equation (5). If the rotor precession 
has a constant radius and rate. /is a constant and the motion 
is linear. However, for general motion, / is variable and a 
function of the motion. 

After performing the asymptotic least-squares curve fit to 
Fig. 2, the resulting coefficients of equation (6) are given in 
Table 2. 

In summary , the following equations represent how the over- 
all coefficients for the model are defined. 

M= M r +2(M S + M l ) 

C=2(C,+Q) + C 

c=2(c, + c t ) (g) 

K = K th + 2(2AT, + Kt) + 2tf* = K + 2K* 

k = 2(k,+k L )+2k’=k + 2k\ 


ramcuiar coefficients are doubled because there are two 
seals and two leakage paths in a double suction pump. The 
quantity, C, must be added to the model to yield reasonable 
stability. If the eigenvalues of the overall system are obtained 
without adding C, the system is found to be unstable. There- 
fore, by adding another damping term (C= 10,422 Ns/m) 
reasonable stability, which is consistent with operating pump 
experience, is ensured. The following equation represents the 
complete model used to model the forces on the rotor. 
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Table 3 Coefficients for overall model 


Table 4 Natural frequencies of free vibration for the pump model with 
and without K '• and Jr* 



M 

C 

c 

K 

k 


*9 

Ns/m 

Ns/m 

MN/m 

MN/m 

.5 



2097.1 

2.8827 

1.5549 

.6 

98.298 

18830 

RMI 


1.5494 

.7 

98.361 

16438 

BdEMEB 

2.3363 

1.5184 



Fig. 6 Real part of the eigenvalues for the forward-whirl mode 



Whirl F requ en cy Ratio (NP) 

Fig. 7 Imaginary part of the eigenvalues for the forward-whirl mode 
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U] r c ci ( 


[ 0 M 

w L- cC Ji 
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r K(f) k(f) 

[-k(f) K(f) 




Table 3 contains the numerical values for this analysis. These 
coefficients combine to make the damping for the uncoupled 
system, defined by C/2\J KM, equal to 63 percent, 59 percent, 
and 54 percent of critical damping for the three swirl ratios 
shown, respectively. Again, ft* and £* are illustrated in Fig. 
5 and are defined analytically with cubic-spline curve fits. 


Natural Frequency Determination 

The characteristic equation for Eq. (9) is 

[Ms 1 + Cs + K(f)'? + [cs + k(f )] 1 = Q, (10a) 

where s = s/«,. For each whirl-frequency ratio /, two pairs of 
complex conjugate roots of the form 

I-a/iyX,; /= 1,2 (10b) 

can be calculated corresponding to forward and backward 
whirling modes. The concern here is finding the natural fre- 


UtL 


fa9EE1113il3 

E : - 

Percent 

Decrease 

.5 

.747 

.727 

2.68% 

.6 

.777 

.724 

6.82% 

.7 

■II 

.715 



Tabla 5 Real parts of the eigenvalues for the model with and without 
K • and k m 



Real part with 
and k * 
1/sec 

Real part without 
K * and k* 

1 /sec 

Percent 

Decrease 

.5 

-55.9 

-60.3 

7.87% 

.6 

-31.1 

-47.2 

51.80% 

.7 

-13.8 

-34.5 

150.00% 


quency for forward whirling motion. Figure 7 illustrates the 
eigenvalue solutions to Eq. (10a) versus/. 

Free whirl of the rotor occurs only when the imaginary part 
of an eigenvalue equals the whirl frequency ratio /in Fig. 7; 
i.e., X 2 =/. The three points at which the X 2 =/line intersects 
the X(/) curves defines the forward-whirl natural frequencies 
for the three inlet swirl ratios. The real part of the eigenvalue 
at this natural frequency defines the rotor’s damping. 

To find the “natural frequency” and damping of free whirl, 
the following iterative technique has been used. First, a whirl 
frequency ratio is assumed. From this value, the corresponding 
values of K and k are obtained from the data shown in Fig. 
5 from which the eigenvalues are found using Eq. (10). These 
first two steps are equivalent to finding a point on the curves 
of Figs. 6 and 7. The imaginary part of the forward-whirl 
eigenvalue is compared to the whirl frequency used. If they 
are the same, the natural frequency at which free vibration 
takes place is defined. If they are different, the imaginary pan 
of the forward-whirl root becomes the assumed whirl frequency 
ratio, and the same steps are carried out until convergence 
occurs. Note that this calculation procedure yields the damped 
natural frequency for the system. 

Tables 4 and 5 show results of the above algorithm. Also 
included are the eigenvalues of the same model neglecting K* 
and k *, which correspond to a pure quadratic approximation 
of the leakage path forces as in equations (4). The eigenvalues 
illustrate at least two important conclusions about the effects 
of the values of K* and k *. First, the natural frequency in- 
creases when K* and k * are included as well as when the swirl 
ratio is increased. The variation of K* in Fig. 5 explains these 
observations. More importantly, the percent differences be- 
tween the two models show that the values of K* and k* have 
very little effect on natural frequency. Second, stability is re- 
duced when A'* and k * are included and when the swirl ratio 
is increased. The variation of k* in Fig. 5 is the cause of these 
observations. Finally, the percent differences between the two 
models shows that K * and k * can cause appreciable reductions 
in relative stability. 

To show that the stability and frequency of free vibration 
are determined only by the root obtained from the iterative 
approach described earlier, the equations of motion (9) were 
integrated in a series of transient, free-vibration simulations. 
Initial conditions were an initial displacement of . 127 mm (.005 
in) and velocity equal to Rlw s (synchronous whirl). Figure 8 
includes the resulting orbit and a time history of the whirl 
frequency ratio for free vibration and an inlet swirl ratio of 
0.7. As shown, the rotor executed a well-damped spiral orbit. 
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Fig. 8 (a) Rotor orbit during free vibration 



Fig. 8(b) Tima history of rotor whirl frequency ratio during free vibration 



W>url Frequency Raw (ND> 

Fig. 9 Real part of the eigenvalues for the forward-whirl mode with a 
12.5 percent reduction in direct damping 


The whirl frequency ratio of the orbit was exactly 0.8 T which 
corresponds to the natural frequency found earlier. The other 
two swirl ratios yield similar results. 

To verify that the stability of the model is described only 
by the real part of the eigenvalues obtained with the iterative 
approach described above, the direct damping coefficient was 
reduced by 12.5 percent producing a range of unstable eigen- 
values as shown in Fig. 9 where the real part of the forward- 
whirl root is positive between whirl frequency ratios of 0.41 
and 0.57. The natural frequency of this altered model corre- 
sponds to a stable eigenvalue and is .825 w,. Integration of the 
equations of motion were altered to include a harmonic forcing 
function with an amplitude of 5N which forced the rotor in a 
circular direction at a frequency of .5w s . Figure 10 includes 
the resulting orbit, whirl frequency ratio time history, and a 
Fast Fourier Transform (FFT) of the jc-coordinate, respec- 
tively. Even though the forcing frequency was within the un- 
stable range of frequencies, the figures show that the system 
was stable in that its orbit did not grow without bound. The 
loops and limit cycle behavior are the results of the nonlin- 
earities in the model. The FFT shows that the steady state orbit 
consisted of response at both the forcing frequency, .5w„ and 
the natural frequency, .825 w,. Other forcing frequencies were 


0.02 


o 

£ 



- 0 . 02-1 1 1 1 * 1 

- 0.02 - 0.01 0,00 0.01 0.02 
x-CocMTiinaie/ShreMid Inlet Clearance (ND) 

Fig. 10(a) Rotor orbit during harmonic excitation at 50 percent of shaft 
speed and 12.5 percent reduction In direct damping 



Fig. 10(b) Time history of rotor whirl frequency ratio during harmonic 
excitation at 50 percent of shaft speed and a 12.5 percent reduction in 
direct damping 



Fig. 10(c) FFT of the response of the rotor during harmonic excitation 
at 50 percent of shaft speed and a 12.5 percent reduction in direct 
damping 


also tried. Using a stable value of 0.7w„ for example, the 
response did not behave in the same mann er; instead, the orbit 
reached a steady state circular path at the excitation frequency 
as would be expected from a linear system. However, forcing 
at .375w„ the characteristics of the results from forcing at .5w s 
occurred again, as can be seen in Fig. 11, which again contain 
the orbit, time history of the whirl frequency ratio and an FFT 
of the response. Only forcing frequencies below and within 
the unstable zone resulted in this behavior. Only when the 
direct damping coefficient was reduced until the natural fre- 
quency corresponded to an unstable eigenvalue did the re- 
sponse become unstable in a linear sense, where the orbit grew 
without bound at the natural frequency. The fundamental 
result from this analysis is a verification that a range of whirl 
frequencies yielding eigenvalues with positive real parts does 
not cause the system response to grow without bound unless 
it includes the system natural frequency. In addition, the non- 
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i-Coonhnaie^hroud Inlet Geennce (ND) 

Fig. 1 1(*) Rotor orbit during harmonic excitation at 37.5 percent of shaft 
speed and a 12.5 percent reduction In direct damping 



Fig. 12 Frequency response: whirl amplitude 



Fig. 11(b) Time history of rotor whirl frequency ratio during harmonic 
excitation at 37.5 percent ol shaft speed and a 12.5 percent reduction 
In direct damping 



Whirl Frequency Rjiuo (ND) 



Fig. 11(c) FFT of the response of the rotor during harmonic excitation 
at 37.5 percent of shaft speed and a 12.5 percent reduction In direct 
damping 


linearities of the system can cause a limit cycle to occur under 
certain conditions. 

Finally, the critical speed of the rotor cannot be calculated 
unless force curves are available at other shaft speeds. In this 
case, the natural frequency is calculated with the iterative 
method just described, but at a number of shaft speeds. A 
plot of natural frequency as a function of shaft speed is then 
drawn, analogous to “critical speed maps” in fluid-film bear- 
ing analysis. The critical speed is defined as the speed where 
the natural frequency is equal to the shaft speed. 


Frequency Response 

As discussed earlier, the following are the equations of mo- 
tion of the model 


•(- 


K(f) k(f) 
k(f) K(f) 


8 - 


.fcosfl/) 
[sinQ/J ’ 


(ID 


Fig. 13 Frequency response: phase angle 


where the right-hand side represents a harmonic forcing func- 
tion appropriate to frequency response analyses. As in tradi- 
tional methods, the transient solution is assumed to decay with 
time, and only the steady-state response of the rotor is con- 
sidered. Obviously, frequency response solutions are only valid 
when the system is stable, precluding the type of response 
illustrated in Figs. 10 and 1 1 . The forcing frequency is defined 
by Q which has been used to represent whirl speed. Since the 
steady-state motion is assumed to occur at the same frequency 
as the excitation, this should not cause confusion. Note that 
the model is now linear because the whirl frequency ratio and, 
hence, the stiffnesses are constant. Multiplying the second 
equation by the complex number /, the equations can be com- 
bined into one in terms of the complex variable z. 

Mz+ (C-ic)z+ (K-ik)z = Fe ia ‘. (12) 

Assuming steady-state motion of the following form 

z = Re/*, (13) 

where R denotes the whirl amplitude, and <t> is the phase angle 
yields 


J? 

<f> = tan “ 1 


K 

-A*tf + cG] 2 +[CQ--A :] 2 
f Ar-CQ ] 
K-MQ 2 + cQ * 


(14) 


Equations (14) were solved using the present model data, 
and the results for a swirl ratio of 0.7 are shown in Figs. 12 
and 13. Also included are the same results for the model with 
K * and k* neglected. Again, in obtaining frequency curves as 
shown in Figs. 12 and 13, the stiffness terms must be calculated 
at the corresponding whirl frequency ratio. 

The swirl ratio of 0.7, used in Figs. 12 and 13, exhibited the 
largest variation in amplitude and phase at the whirl frequency 
ratio of .80, the system natural frequency. An amplitude dif- 


r* + m 






fercncc of about 60 percent is shown between the curves, sug- 
gesting that K * and k * have significant influence on the system 
primarily at this frequency. The swirl ratios of 0.6 and 0.5 
exhibited similar but decreasing results in both the amplitude 
and phase curves, having differences in amplitude at the peak 
value of 35 percent and 6 percent, respectively. In fact, the 
curves were almost identical in the case corresponding to a 
swirl ratio of 0.5. Note that these resonance peaks occurred 
in all cases at the corresponding natural frequencies found in 
the previous section, further verifying the preceding results. 
The difference in the response to synchronous imbalance ex- 
citation (/= 1) between the two models is small, and this sug- 
gests that the leakage path forces have negligible effects on 
imbalance response. However, the presence of a subsynchron- 
ous excitation near the natural frequency can result in large 
increases in response amplitude when k* is included. 

Summary and Conclusions 

The impeller-shroud forces have been separated into (a) con- 
stant (whirl frequency-independent) stiffness, damping and 
mass coefficients and ( b ) direct and cross-coupled added stiff- 
ness coefficients which are functions of whirl frequency. With 
this approach, the model can be analyzed using traditional 
techniques with a few modifications. A new iterative technique 
was used in the determination of the rotor’s natural frequency 
and damping. The whirl-frequency-dependent stiffness terms 
were seen to be convenient and effective models for the im- 
peller-shroud force nonlinearities. 

When the values of the added stiffness coefficients, K and 
k* t were included in addition to the frequency-independent 
coefficients, a number of effects resulted. First, the natural 
frequency was seen to increase by a small amount, suggesting 
that AT* could be neglected when only the natural frequency 
is important. In addition, the stability is seen to decrease when 
k* is included. This effect is appreciable, so stability calcu- 
lations should include k\ The frequency response analysis 
showed synchronous response to imbalance is increased by the 
stiffnesses although by only a small amount. Appreciable dif- 
ferences in response due to subsynchronous excitation occur 
only at or near the natural frequency and become negligible 
at low values of the inlet swirl ratio. 

The nonlinearity introduced by the whirl-frequency-de- 


pendency of the force coefficients resulted in limit cycles and 
nonsynchronous response in some cases. 


References 


Adtkins, D., 1976, “Analyses of Hydrodynamic Forces on Centrifugal Pump 
Impellers ” M.S. Thesis, California Institute of Technology, Pasadena. 

Bolletcr, U., Wyss, A.. Welte, I., and Sturchler, R., 1985, “Measurements 
of Hydrodynamic Interaction Matrices of Boiler Feed Pump Impellers,” ASME, 

Boileter, U., Leibundgut, E., Sturchler, R., and McCloskey, T., 1989, ^ Hy- 
draulic Interaction and Excitation Forces of High Head Pump ImpeUers,” pre- 
sented at the ASCE/ASME Pumping Machinery Symposium, July 9-12, 1989. 

Childs, D., Dressman, J., and Childs, S., 1980, “Testing of Turbulent Seals 
for Rotordynamic Coefficients," Proc. of Workshop on Rotordynamic Insta- 
bility Problems in High Performance Turbomachinery, Texas A&M University. 
NASA CP-2133, pp. 121-38. 

Childs D., 1981, “Convergent-Tapered Annular Seals: Analysis for Rotor- 
dynamic ’Coefficients,” Symposium Vol.. Fluid/Structure Interactions Turbo- 
mach., ASMW Winter Annual Meeting, pp. 35-44. 

Childs, D., 1982a, ‘'Dynamic Analysis of Turbulent Annular Seals Based on 
Hirs’ Lubrication Equations,” ASME Paper No. 82-Lub-41. 

Childs, D., 1982b, “Finite-Length Solutions for Rotordynamic Coefficients 
of Turbulent Annular Seals,” ASME Paper No. 82-Lub-42. 

Childs, D.. and Dressman, J.. 1982c, “Testing of Turbulent Seals for Ro- 
lOf dynamic Coefficients, “ Proceedings of Workshop on Rotordynamic Insta- 
bility Problems in High Performance Turbomachinery , Texas A&M University, 

NASA CP-2250, pp. 157-71. , u unrvrp 

Childs, D. t and Moyer, D., 1985, “Vibration Characteristics of the HPOTP 
(High Pressure Oxygen Turbopump) of the SSME (Space Shuttle Main Engine),” 
ASME Journal of Engineering for Gas Turbines and Power, Vol. 107, pp. 152- 
159. 

Childs, D., 1987, “Fluid-Structure Interaction Forces At Pump-Impeller- 
Shroud Surfaces for Rotordynamic Calculations,” Texas A&M University, Me- 
chanical Engineering Dept., ASME Vibration Conference, Boston. Mass. 

Childs, D., 1990, “Centrifugal-Acceleration Modes for Incompressible Fluid 
in the Leakage Analysis Between a Shrouded Pump Impeller and Its Housing,” 
ASME Journal of Vtoatton and Acoustics, May 1990. 

Franz R and Arndt, N., 1986, “Measurement of Hydrodynamic Forces on 
the impeller of the SSME.” Report No. E249.2, California Institute of Tech- 
nology, Pasadena. . w . . . 

Jcffcott H 1919, “The Lateral Vibration of Loaded Shafts in the Neigh- 
borhood of a Whirling Speed-The Effect of Want of Balance,” The Philo- 
sophical Magazine, Vol. 6, No. 37, pp. 304-314. 

Jery, B.. Acosta, J., Brennen, C. ( and Caughey, T., 1984, Hydrodynamic 
Impeller Stiffness, Damping and Inertia In the Rotordynamics of Centrifugal 
Flow Pumps,” Proceedings of Workshop on Rotordynamic Instability Problems 

u:~l Turhnmnrhinsrv . Texas A&M University, NASA CP-2338, 


Massey, I. C., 1985, “Subsynchronous Vibration Problems in High Speed 
Multistage Centrifugal Pumps,” Proceedings of the 14th Turbomachinery Sym- 
posium , Turbomachinery Laboratories, Texas A&M University, pp. 11-16. 


Journal of Vibration and Acoustics 


OCTOBER 1991, Vol. 113/515 


i 

; 

I 


( 1992b) 




Reprinted from March 1992, Vol. 114, Journal of Fluids Engineering 




D. W. Childs 

Turbomachinery Laboratories, 
Mechanical Engineering Department, 
Texas A&M University, 
College Station, Texas 77843 


Pressure Oscillation in the 
Leakage Annulus Between a 
Shrouded Impeller and Its 
Housing Due to Impeller- 
Discharge-Pressure Disturbances 

An analysis is presented for the perturbed flow in the leakage path between a 
shrouded-pump impeller and its housing caused by oscillations in the impeller- 
discharge pressure. A bulk-flow model is used for the analysis consisting of the path- 
momentum, circumferential-momentum , and continuing equations. Shear stress at 
the impeller and housing surfaces are modeled according to Hirs ' turbulent lubri- 
cation model . In the present analysis, perturbations of the impeller discharge pressure 

are used to excite the fluid annulus. The circumferential variation of the discharge 
pressure is expanded in a Fourier series up to order n u where n x is the number of 
impeller blades. A precession of the impeller wave pattern in the same direction or 
opposite to pump rotation is then assumed to completely define the disturbance 
excitation. Predictions show that the first (lowest-frequency) "centrifugal-accel- 
eration ” mode of the fluid within the annulus has its peak pressure amplitude near 
the wearing-ring seal . Pressure oscillations from the impeller can either be attenuated 
or ( sharply ) magnified depending on: (a) the tangential velocity ratio of the fluid 
entering the seal, (b) the order of the fourier coefficient, and (c) the closeness of 
the precessional frequency of the rotating pressure field to the first natural frequency 
of the fluid annulus, and (d) the clearance in the wearing-ring seal. 


Introduction 

The present work is stimulated by experiences with the SSME 
HPFTP (Space Shuttle Main Engine, High Pressure Fuel Tur- 
bopump) wearing-ring seals. A stepped, 3-cavity, tooth-on- 
rotor, labyrinth-seal design is used. The stator for the seal is 
made from KEL-F, a plastic that is somewhat similar to nylon. 
In some cases, post-test inspection of the stator element has 
revealed that interior points in the stator material have melted 
and then resolidified, despite being in contact with liquid hy- 
drogen. One hypothesis for this exceptional outcome was that 
the material had been subjected to cyclical stresses which gen- 
erated heat due to hysteresis. Because of poor conduction 
properties of the material, the heat could not be dissipated, 
the temperature rose, and melting resulted. “What pressure 
oscillations are driving the cyclical stresses?”, is an obvious 
question in reviewing this scenario. The present analysis ex- 
amines “centrifugal-acceleration” modes, arising between the 
impeller shroud and its housing and driven by pressure oscil- 
lations from the pump, as an answer to this question. 


'The work reported herein was supported by NASA Marshall Space Flight 
Center under contract NAS 8-37821; contract technical monitor: James Cannon. 

Contributed by the Fluids Engineering Division for publication in the Journal 
or Fluids Engineering. Manuscript received by the Fluids Engineering Division 
August 20, 1990. 
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Figure 1 illustrates an impeller stage of a multi-stage cen- 
trifugal pump. Leakage along the front side of the impeller, 
from impeller discharge to inlet, is restricted by a wearing-ring 
seal, while leakage along the back side is restricted by either 
an interstage seal or a balance-piston discharge seal. The pres- 
ent analysis considers perturbed flow in the leakage paths be- 
tween the impeller-shroud surface and its housing. 

Prior analyses by the author of those annulli have been 
concerned with lateral (1987, 1989) and axial (1990a) reaction 
forces developed by the impeller shrouds as a consequence of 
harmonic clearance changes due to impeller motion. These 
analyses have been based on “bulk- flow” models which neglect 
the variation in the dependent variables across the fluid film. 
The model consists of the path and circumferential momentum 
equations and the continuity equations. 

The analyses cited have yielded force and moment coeffi- 
cients due to impeller motion but have also predicted “reso- 
nance” phenomena, which are caused by the centrifugal- 
acceleration body forces present in the path momentum equa- 
tions. An algorithm was developed (1990b) to calculate the 
complex eigenvalues and eigenvectors associated with these 
resonances. In the present analysis, the harmonic response of 
the flow within the annulus is examined due to time and cir- 
cumferential variations in the discharge pressure of the im- 
peller. 
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Geometry and Kinematics 

Figure 1 illustrates the annular leakage paths along the front 
and back sides of a typical shrouded impeller of a multistage 
centrifugal pump. The present discussion concentrates on the 
flow and pressure fields within the forward annulus; however, 
the analysis can also be applied to the rear annulus. 2 As il- 
lustrated in Fig, 2, the outer surface of the impeller is a surface 
of revolution formed by rotating the curve R = R(Z) about 
the Z axis. A point on the surface may be located by the 
coordinates Z, R (Z) , 0. The length along the curve R (Z) from 


Y 



the initial point R h Z, to an arbitrary point R, Z is denoted 
by S and defined by 



In the equations which follow, the path coordinate S and 
angular coordinate B are used as independent spatial variables. 
The coordinates Z, R defining the impeller surface are ex- 
pressed as parametric functions of S, i.e., Z(S), R(S). The 
length of the leakage path along the impeller face is defined 
by 



Although ihe leakage flow is normally up the backside of all impellers except 
the last impeller, the governing equations would continue to be valid irrespective 
of the flow direction. 


Governing Equations 

Returning to Fig. 2, the path coordinate 5 and circumfer- 
ential coordinate RB are used to locate a fluid differential 


Nomenclature 


A 2s > A} St ^26, Ai0 = 
b = Vf/Ru = 

C* = 

Q = 
Cr = 

/= 0/a> = 
/• - 


r 1 


”1 

(ni-n 2 ) 


h - H/C t = 
H = 

L s = 

n - 

P=p/ P vj = 

A(iV) = 

p t (6,t) = 


coefficients introduced in Eq. (12) 
nondimensional velocity ratio 
discharge coefficient for the exit wear- 
ring seal introduced in Eq. (7) 
initial (5 = 0) clearance (Z,) 

exit seal clearance (L) 

nondimensional precession frequency 
nondimensionalized precession fre- 
quency yielding a maximum response 
pressure 

positive, dominant, nondimensionaI , 
precession frequency predicted for an 
impeller with n x blades in a diffuser 
with n 2 blades 

nondimensionalized clearance 
clearance between impeller shroud and 
housing (L) 

leakage-path length, defined by Eq. 

( 2 ), (L) 

order ofPourier coefficient, intro- 
duced in Eq. (13). 

nondimensional static fluid pressure 
prescribed annulus supply pressure 
(impeller exit pressure) 
prescribed annulus exit pressure 
(impeller inlet pressure) 


P 

R 


Ri 

r-R/Ri 
R s = 2HU s /v 


S 


s = S/L s 


= fluid static pressure ( F/L 2 ) 

= radial coordinate (L) 

= initial (s = 0) radius ( L ) 

= nondimensionalized radial coordinate 
= path-velocity Reynolds number 
= path coordinate introduced in Eq. (1), 

= nondimensionairzed path length 


T~ L S /Vj = representat ive transit time for flui d 
traversing the leakage path (7) ”*■ : 

« 5 = U s /Vi = nondimensionalized path fluid velocity 
u e = U e / Rju = nondimensionalized circumferential 
fluid velocity 

Vi = initial (5 = 0) path fluid velocity 
€ = perturbation coefficient 
o) - pump running speed (T -1 ) 
ft = excitation frequency (T~ *) 
p = fluid density (M/U) 

B = circumferential coordinate 
a s = normalized friction factors, defined by 
Eq. (11) 


7 = rf = nondimensionalized time 

£ = entrance-loss coefficient introduced in 
Eq. (7) 

v = kinematic viscosity (L 2 T ~ ') 

Overbars denote complex variables; see Eq. (16). Subscripts 0 
and 1 denote zeroth- and first-order solutions, respectively. 
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element of thickness H(S, $, t). From Childs (1987), the con- 
tinuity equation can be stated 


dH £ 
dt + dS 


(UsH)+ ^Ve iU<>H) + 


H\ dR „ n 
/?) dS Us 0 


(3 a) 


where U s and U e are the path and circumferential bulk-velocity 
components, respectively. The path and circumferential mo- 
mentum equations are stated 


„U}dR 
H dS~ Ta + Ts ' pH R dS 


H dP 

~r de~ T9!+T9r 


+ pH 


dUs dUs Ue 3(Js 
dt + dd R + dS 



m 


+ pH 


dU„ dUg Ue dU e 
dt + d6 R + dS 


U s + 


U»Us BR\ 
R dSj 


(3c) 


Following Hirs’ approach (1973), the wall shear-stress defi- 
nitions in these equations can be stated 


T a = - pU)RT[\ + (U e /U s ) 


mj+ I 




nr i mr * 1 

T S r = - pU)R?'\ 1 + [(Ue-R^/U,) 1 ] — 

ns ms * 1 

T*s = y pU e U s RT[ 1 + (U 9 /U s ) 2 ] 2 

nr mr + 1 

T<,r=jpU s (U e - Rui)RT r {l + [(U e -Ru>) /U s } 2 1 2 (4) 

where 

R s = 2H\JJv (5) 

Nondimensionalization and Perturbation Analysis 

The governing equations define the bulk-flow velocity com- 
ponents ( U s , U e ) and the pressure Pas a function of the co- 
ordinates ( R8 , S) and time, t. They are conveniently 
nondimensionalized by introducing the following variables 

u s= U s /V„ U) = U e /Rju>, p = P/p Vf 
s = S/L s , r= R/R„ b = Vj/Rjtp 

T=oit, T= L,/Vi (6) 

The present analysis examines the changes in (u s , u e , p) due 
to changes in the impeller’s discharge or inlet pressure. Fol- 
lowing conventional notation, pressure drops at the annulus 
inlet and exit are stated 

PAd.t) - P(0,d,t) = p( 1 + l)U)(0,6,t)/2 

P(L S , 6, t)~PA6,t) =pC dt U](L s ,6,t)/ 2 (7) 

Note specifically that the (upstream) supply and (downstream) 
exit pressure are now functions of time. Assume that the os- 
cillations consist of a small perturbation of the form 

P s (0,t) = P„> + (P st (6,t) ,P e (6,t) = P# + eP tl (6,t) 

and introduce nondimensional variables to yield the following 
zeroth 

A>(0)=Ao-( l-|)/2 

P(V=P a +C de ul 0 (l)/2 (8) 

and first-order equations 

Psi (0>t) -Pi (0,0,0 = ( 1 + £)w„ (0,0, t) 

p x (1,0,0 -P' X (0,/) = C^(\)u sX ( 1 ,0,/) (9) 


The perturbed supply and exit pressures p s ,(0,f), p el (6 y t) can 
now be specified functions of time and provide excitation for 
the perturbed flowfield within the annulus. 

Expansion of the dependent variables of Eq. (3) in pertur- 
bation equations yields: 

Zeroth-Order Equations 


(a) Path-Momentum Equation 


dPo_ \ 

(dr\ 

Ulj 

(o r + 

_±dho I dr 


ds r 

\ds) 

* l 

( ^ ) 

h 0 ds r ds 

|«4=o (io<7) 


(b) Circumferential-Momentum Equation 


du eo u 9o dr 

*57 + — ^ + l°r ( u 0o - r) + CsUeo] /2 = 0 (10 b) 

(c) Continuity Equation 

rh 0 u so = 1 (10c) 

The quantities a s and a r are defined by 

a, = (L s /H 0 )\ st o r = (L s /H 0 )\ r (1 1) 

where \ s and X r are the dimensionless stator and rotor friction 
factors: 

ms+ I 

K = nsR™{ 1 + (u ta /bu so ) 2 ] 2 

mr + | 

A, = nrRZ { 1 + l(u 6o - r) / bu so ] 2 ) ~ 

The continuity equation has been used to eliminate du so /ds 
from Eq. (10a). The momentum equations define the pressure 
and velocity distributions for a centered impeller position. They 
are coupled and nonlinear and must be solved iteratively. The 
initial condition for u $0 ( 0) is obtained from the exit flow con- 
dition of the impeller. Zeroth-order pressure boundary con- 
ditions are provided by Eq. (8). 

Figure 3 illustrates the pump-impeller and shroud geometry 
used by Bolleter et al. (1987) in their test program for radial 
force coefficients. Their tests were at best efficiency point 
(BEP) with the pump running at 2000 rpm, while developing 
68m of head and 1301/s of flow rate. The impeller has seven 
blades and an impeller exit angle of 22.5*. The test fluid is 
water at 26.6 *C. For the present study, A P across the impeller 
is assumed to be 70 percent of the total head rise of the stage. 
Based on pitot-tube measurements, impeller-exit-tangential ve- 
locity is about 50 percent of the impeller discharge surface 
velocity; hence, u m 3 0.5. 
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Table 1 Zeroth-order-solution results; C, * 3.5mm, C r 
0.36mm _ 


u eo ( 0 ) 

RfO 

C dt 

m (kg/sec) 


0.5 

9377 

2.068 

4.448 


0.6 

8907 

2.098 

4.225 


0.7 

8426 

2.130 

3.997 


«5i = w 5 | C +yw 5 | 5 , Ue\ = u eic+J u »\s P\-P\c+JP\s 0 6 ) 

reduces these real equations to three, complex, ordinary dif- 
ferential equations 


Table 1 provides the zeroth-order solutions for new (original- 
clearance) exit wearing-ring seals. The R& values on the order where 
of 10 000 are low in comparison to the circumferential Reyn- 
olds number Reo = THV^p which varies along the path, but 
is on the order of 250,000. These leakage results were obtained 
iteratively starting with guessed values for the seal inlet values ^ = 
of pressure and tangential velocities and yielding an mitia 
estimate for C de . With this estimate for C* t the leakage is 
calculated through the annulus which yields new inlet condi- 
tions for the seal. Solutions are “bounced” back and forth 
between the seal and the annulus until the same leakage value 
is obtained for both flow paths (four-place accuracy). 


£ 

ds 


B 

A^ u so 

A i5 - BUsq+jTT 



1 

M 


+ [A ( n,/,s )] < 
1 1 

L P'J 


(17) 


- jnuit/r 


( A 26 4* T) / UsQ 

l nuT \ 


-J 


0 

nb 

rwso 

0 


1 dr 

* rds h 0 ds' \ r 


(18) 

(19) 


First-Order Equations 

(a) Path-Momentum Equation 

Uq\A2s + W 5 |j435 


Since there is no right-hand side to Eq. (17). the homogeneous 
solution is the complete solution and can be stated as follows 
in terms of the transition matrix and initial conditions 


r t 3uji . t Ue ° d — i u 

+ [“ r 'a7 r r de +u 

(b) Circumferential-Momentum Equation 


dUS { 

ds 


= 0 (12a) 


u$\ 

UJ 


= l4>(n,/,5)] 


/ M J l(0)'j 

«.i(0)j 

fp,(0)J 


( 20 ) 


b-jr- ~TT + 8 + U s\A 

Rj r off 


+ \ wT ^ + uT ^ + ujM=0 (12 b) 

3t r au os 

(c) Continuity Equation 

du^ wTdus I A4 (12c) 

ds r dd 11 \r ds h a ds j 

New coefficients in these equations are defined in, Childs 

(1987). 

Solution Procedure: First-Order Equations 

The functions p sl (6, /), p„(0, 0 provide the boundary ex- 
citation for the first-order equations. The general form for the 
excitation takes the form 

p sl (9,t) =e ja, (Psic cos nd + p ils sin nd) 

P'A0't)=e Ja '(P'U cos nd+P'is sin nd) (13) 

where n can reasonably be expected to vary from zero (plane 
wave) upwards through multiples of the number of blades in 
the impeller. The form ofEq. (13) suggests that the 6 variation 
in boundary pressures is defined in an impeller-fixed coordinate 
system, which is precessing at the frequency fi. .... 

The 6 and time dependency of the dependent variables is 
eliminated by assuming the comparable, separation-of-van- 
able, solution format 

u sl = e Jfr (u, u cos nd + u sU sin nd) 

Ue\=e jfr (u tu cos nd + u « u sin nd) 

Pi=e J/T (Pu cos nd+p ls sin n8) (I 4 ) 

where the coefficients are solely functions of s, and 

is the normalized precession frequency. Substituting into Eqs. 
(12) and equating like coefficients of cos nd and sin nd yields 
six first-order equations ins. Introducing the complex variables 


The inlet initial condition «,,(0) is set equal to zero, and 
calculation of 1^,(0) and p t ( 0) in terms of the specified bound- 
ary conditions is the immediate problem at hand. Substitution 
from Eqs. (13) and (14) into Eq. (9) yields 

Ps\ — Pi(0) = (1 + «)«,.(<» 


Pl(l) ~ Pel — ^-rfpUjo(l)Uji(l) 

where 

Pa =Ps\c+jPsis 

Pel = Pelc jPels 

From Eq. (20) 

« J i(l) = *ii(1)«ji( 0) + *u(l)Pi(0) 

Pl(l) = 4>3|(1 )«j|( 0) + ^3j(l)Pl(0) 

Hence, from Eq. (21) one obtains 

[ (l+O 1 

|_^3l(l — Cte«jo(l)^ll(l) 

x . fa.) 

(R(0)j {Pel) 


( 21 ) 


( 22 ) 


(23) 


(24) 


fi J iviuj 

r Zii Zu 

|_ Z21 Z22 {pel) 


(25) 


Inversion of this equation yields 

(« 1 i(0)') 

[pi(0)j 

p s[ and Pti cannot be specified independently, and a relation- 
ship between the two cannot be established without a knowl- 
edge of the fluid system beyond the current terminating orifices. 
For the purposes of this discussion, the arbitrary choice 

fci-1. P< i =0 

is made to examine the influence of pressure perturbations at 
the impeller exit (annulus inlet). The resulting set of initial 
conditions for Eq. (20) is then 


rc^on rz,r) 
U(0)j (z 21 j 


u*i(0) = 0 


(26) 


The complete solution along the impeller is found by evaluating 
Eq. (20) for s€[0, 1J. 
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J uqo(O) = 0.5 



S 

Fig. 4 p,(s) complex eigenvector from Childs (1990b); n = 1 


Numerical Results. Childs’ ( 1 990b) analysis yields complex 
eigenvalues a and eigenvectors u s] (s ), w*i(s), Pi(s) for the 
system modeled by Eq. (12). Figure 4 illustrates the amplitude 
and phase for the First (lowest-natural-frequency) eigenvalue 
for n = 1; w w ( 0) = 0.5, 0.6, and 0.7. Observe that the peak- 
pressure amplitudes lie near the wearing-ring seal, which is 
consistent with the internally-melted HPFTP wearing-ring-seal 
results cited earlier. Forced, harmonic-response solutions de- 
veloped in this study due to impeller-discharge-pressure oscil- 
lations also show the largest pressure oscillations to occur near 
the wearing-ring seal (s = 0.95); hence, results presented here 
focus on Pi(0.95). This peak-pressure-oscillation location is 
very near the exit-wearing-ring seal of Fig. 3. The first question 
to be addressed here is, “How does pi(0.95) depend on n,/, 

and C,?” _ ^ 

Figure 5 illustrates the amplitude and phase of p^O.95) versus 
/for C, = 0.36mm, n = 5, and u w = 0.5. As will be explained 
later, the choice n = 5, arises because of the number of impeller 
and diffuser blades used in Bolleter et al.’s (1987) pump. The 
phase results indicate that numerous resonances exist for po- 
sitive values of /. However, only the first resonance experiences 
significant amplification. The remaining fluid modes are heav- 
ily damped. 

Figure 6 illustrates lp,(0.95)l for C, = 0.36mm, n = 5, and 
„ w (0) = 0.5, 0.6, and 0.7. The peak-response frequency in- 
creases as u«o(0) is increased from , 0,5 to 0.6, and a secondary 
peak appears around / = 2.5. Increasing u M ( 0) from 0.6 to 
0.7 causes an additional peak to appear. 

Figure 7 illustrates lpi(0.95)l for C, = 0.36mm, ^so(O) = 
0.5, and n = 0, 1, 3, 5, and 7. The response is heavily damped 
for n = 0, rises sharply as n is increased to one, but then 
remains relatively constant as n ranges upwards over 3, 5, and 
7. Figure 8 repeats the results of Fig. 7, except for worn clear- 
ances; i.e., C, = 0.72mm. Comparisons of Figs. 7 and 8 show 
that doubling the clearances reduces pressure amplification and 



NOND I MENS I ON AL EXCITATION FREQUENCY 

Fig. 5(a) ip, (0.95) I versus f; C, « 0.36mm, (W0) = 0.5, n = 5 



-4 . 0 - A. 'O V. W o.i0 i kj.k 

nondimensional excitation frequency 


Fig. 5 (b) Phase £,(0.95) versus /; C r * 0.36mm, u^O) = 0.5, n «= 5 



and 0.7 


slightly elevates the peak-amplitude-excitation frequency/*. 
Table 2 shows /* versus n for the new(C r = 0.36mm) and 
worn sea!s(C r = 0.72mm). Note that/* increases more-or-less 
linearly with increasing n. 

The questions which now arise are: In a real pump, what 
value of n is likely to arise in impeller-pressure-discharge pat- 
terns, and what precession frequency is most likely to be present 
and dominant? Answers to these questions have been provided 
by Bolleter (1988), who presents an analysis for the pressure 
waves developed by the interaction of impeller and diffuser 
vanes or impeller vanes and volutes. For an impeller with n j 
vanes and a diffuser with vanes, Bolleter shows that a ro- 
tating pressure wave is developed around the impeller exit with 
n = |/i t - n 2 l diametral nodes. If n x > n 2 the pressure wave 
rotates in the direction of the pump with the frequency >/ 


Journal of Fluids Engineering 


MARCH 1992, Vol. 114/65 



Table 2 / * (peak-excitation-amplitude frequency) versus n Table 3 Dominant normalized precession frequency f + and peak-ex- 
for w w (0) = 0.5 citation-amplitude frequency / • versus n, and n 2 ; u«(0) = 0.5 


n 

1 

2 

3 

4 

5 

6 

7 

C T - 0.36mm 

0.3 

0.8 

1.4 

2.1 

2.9 

3.6 

4.2 

C, = 0.72mm 

0.3 

0.9 

1.6 

2.4 

3.1 

3.8 

4.4 
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Fig. 7 I p 1 (0.95) i versus f for n = 0, 1, 3, 5, 7; C, = 0.36mm f and u*(0) 
= 0.5 
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Fig. 6 Ifr (0.95) I versus f for n = 0 t 1, 3, 5, 7; C, * 0.72mm, and u«(0) 
= 0.5 

rt. If rt, < n 2 the precessional frequency is -n^/n. Note that 
n, the number of diametral nodes, cited by Bolleter is the same 
n used in Eq. (13) for the pressure excitation. Further, in terms 
of Eq. (13), 0 * ±n l o)/\n l - n 2 \ . Tyler and Sofrin (1962) 
earlier developed this same result in analyzing the noise gen- 
erated by the interaction of a rotor and stator in axial com- 
pressors of gas turbines. 

Bolleter et al.’s pump (1987) used 12 diffuser blades, and 
the impeller of Fig. 3 has 7 blades. Hence, from Bolleter (1988), 
n = 17 - 121 - 5,andQ = -7«/5 « - I.4«. If the impeller 
were mounted in a double volute, n = 17-21 * 5, and Q 
~ 7ti)/5 = 1.4u. From Fig. 5 (n = 5, u w (0) = 0.5), ampli- 
fication for / = - 1,4 and 1.4 is 0.23 and 0.5, respectively. 
Hence, pressure disturbances from the impeller would generate 
pressure oscillations about twice as large in a double volute as 
in a 12-vaned diffuser. However, in either case, because the 
predominant frequency is well removed from the peak-am- 
plitude-excitation frequency/* = 2.9, impeller pressure dis- 
turbances would actually be attenuated by the annulus. 

From Bolleter’s equations, and the results of Figs. 5 through 
7, significant amplification of impeller-discharge-pressure var- 
iations will only arise when the number of impeller blades 
exceeds the number of diffuser (or volute) blades, yielding a 
positive normalized precession frequency. 


f* =/i,/h = -n 2 ) (27) 
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(C r *0.5«mm) 

r 

(C r «O.T3mm) 

8 

1 

7 

1.14 

4.2 

4.4 


2 

6 

1.33 
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Fig. 9 ip,(0.95)l versus f ; n = 3, 0) = 0.5; C, = 0.36 and 0,72mm 


Moreover, for significant amplification within the leakage an- 
nulus, /+ must lie near/*, the peak-amplitude excitation fre- 
quency. Table 3 shows the variation of/ + and /* for various 
combinations of n x and n 2 . The case of n x = 8 (eight-bladed 
impeller) and n 2 = 4 (four-bladed diffuser) yields a close prox- 
imity of/* = 2.1, 2.4 to /* = 2.0; however, this is an 
unrealistic combination. For a practical configuration, the 
nondimensional frequencies/ 4 = I.67and/* = 1.6areclosest 
for n x = 5(five-bladed impeller) and n 2 = 2(double-discharge 
volute). Figure 9 illustrates lpi(.95)[ for u^O) = 0.5 and n 
= 3, confirming the predictions of Table 3. An amplification 
by a factor of 2.6 is predicted for new clearances and 3.2 for 
worn clearances. 

Numerical Uncertainty. The numerical uncertainty issue 
for the results presented concerns the numerical integration of 
Eq. (17). The results presented were obtained with a fourth- 
order Runge-Kutta integrator package using 200 integration 
steps for the interval [0, 1], Repeating these calculations with 
400 integration steps yielded the same results to about three 
significant figures. 
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Summary and Conclusions 

An analysis has been developed and results presented for 
the pressure oscillations in the leakage annulus between a 
shrouded pump impeller and its housing. These pressure os- 
cillations are driven by a circumferential variation of the im- 
peller discharge pressure which can precess either in the same 
or opposite to the direction of rotation. The circumferential 
variation can be modeled with a Fourier decomposition with 
each mode having n diametral nodes across the impeller. The 
peak-pressure oscillations within the impeller are predicted to 
occur near the exit wearing-ring seal in association with a 
centrifugal-acceleration-mode response of the fluid within the 
annulus (Childs, 1988, 1990b). The peak-amplitude-excitation 
frequency increases linearly with n. Using Bolleter’s (1988) 
work which provides a dependency of n and the precessional 
frequency on the number of impeller (*,) and diffuser (/i 2 ) 
blades, situations are presented which can yield large ampli- 
fications (or significant attenuation) of impeller discharge var- 
iations. The occurrence and nature of the pressure oscillations 
are shown to depend on: (a) the tangential-velocity ratio of 
the fluid entering the seal, ( b ) the order of the Fourier coef- 
ficient, (c) the closeness of the precessional frequency of the 
rotating pressure field to the first natural frequency of the fluid 
annulus, and (d) the clearance of the wearing-ring seal. 


The present results suggest an explanation for the internal 
melting observed on SSME HPFTP seal parts. However, given 
liquid hydrogen’s significant compressibility, a more complete 
analysis, including fluid compressibility, is in order. 
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